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The dynamics of cold Bose atoms in driven tilted optical lattices is analyzed focusing on de- 
struction of Wannier-Stark localization and the phenomenon of band collapse. It is argued that an 
understanding of the experimental results requires thorough account for interaction effects. These 
are suppression of the ballistic spreading of atoms for resonant driving (a multiple of the driving 
frequency coincides with the Bloch frequency) and unbounded sub-diffusive spreading of atoms for 
off-resonant driving. 



I. INTRODUCTION 

Since the first realization of an atomic Bosc-Einstein 
condensate (BEC) in 1995 much attention is payed to 
the role of atom-atom interactions in different coherent 
phenomena of the single-particle quantum mechanics. In 
particular, restricting ourselves to BECs in optical lat- 
tices [H, papers 0, [J discuss the effect of atom- atom in- 
teractions on interband tunneling, papers 0-0] on quasi- 
momentum Bloch oscillations in tilted optical lattices, 
papers [§-fll| on the phenomenon of Anderson localiza- 
tion in a disordered or quasiperiodic ID lattice. 

In the present work we study the effect of inter-atomic 
interactions on BEC's dynamics in driven tilted optical 
lattices. In recent years driven lattices were intensively 
studied experimentally with respect to the so-called phe- 
nomenon of the Bloch band collapse p^4T3 | . This term 
comes from the theoretical prediction for the width of 
the quasienergy Bloch band, which can take zero values 
at certain values of the driving amplitude [HI, [r| . This 
effect is also present in driven tilted lattices if the driv- 
ing frequency is commensurate to the Bloch frequency, 
defined by the tilt [l3j . In this work we revisit the prob- 
lem of band collapse in driven tilted lattices, focusing on 
the role of interactions. 

The second fundamental problem we address in this 
paper is the interaction-induced destruction of the 
Wannier-Stark localization. It is known that for van- 
ishing inter-atomic interactions the eigenfunctions of an 
atom in a tilted lattice are localized Wannier-Stark states 
and, hence, any (initially) localized wave-packet remains 
localized during time evolution. It was argued recently 
[13, 13 that, for finite interactions, the time evolution of 
the wave-packet may be unbounded. It should be stressed 
from the very beginning that this regime requires weak 
static forces, not easily accessible in a laboratory exper- 
iment. In the opposite limit of a strong static force, the 
wave-packet dynamics is always bounded, as it has al- 
ready confirmed experimentally [19j . A new generation 
of the cited experiment [ljl] allows to enter the regime of 
weak static forces by using an additional periodic driv- 
ing [2(| Hlj]. This opens perspectives for studying the 
interaction-induced destruction of the Wannier-Stark lo- 



calization experimentally. 

Our theoretical framework is the following single- 
particle Hamiltonian 
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where \l) are Wannier states, J the hopping matrix ele- 
ments, d the lattice period, F the magnitude of a static 
field, F u and u the magnitude and frequency of AC field. 
Follow ing the detection scheme of the laboratory experi- 
ments [3, H3 we are interested in the spacial spread- 
ing of an initially localized wave-packet, which we char- 
acterize by the square root of the wave-packet second 
momentum 



a(t)=[Y J l 2 Piit)~x 2 (t) 



1/2 



x{t)=Y, lP ^) ( 2 ) 



(here P/ (t) are the occupation probabilities of the lattice 
sites, $Z;-Pi(t) = 1)- As known, for vanishing interac- 
tions the system |1| allows an analytical treatment and 
the quantity a(t) can be calculated exactly for an arbi- 
trary initial state [22], [23[. In what follows we borrow 
the relevant equations for the first and the second wave- 
packet momentum from [22| . where the model ((IJ was 
treated by means of dynamical Lie algebras. 

The structure of the paper is as follows. Section |TT] 
is devoted to the dynamics of non-interacting and inter- 
acting atoms in a stationary lattice. We recall essen- 
tials of the unbounded regime and accomplish studies of 
Ref. [13, HU by analyzing the rate of the wave-packet 
spreading in dependence on the strength of atom-atom 
interactions. The case of driven lattices is considered 
in Sec. IIIIl It is shown that the problem of BEC dy- 
namics in a driven lattice can be mapped to that in the 
stationary lattice with properly renormalized static field 
magnitude and hopping matrix elements. This allows an 
understanding of the main features of BEC dynamics in 
driven lattices by referring to the static case of Sec. [TTJ 
The main results of the work are summarized in the con- 
cluding Sec. [TV] 
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II. STATIONARY LATTICES 
A. Single-atom dynamics 

Without driving, F u — 0, the eigenfunctions of the 
Hamiltonian ([1} are localized Wannier-Stark states 

\m) = Y,J l - n (J/dF)\l) , (3) 
i 

(here J v (z) are Bessel functions of the first kind) and 
the spectrum is given by the Wannier-Stark ladder with 
level spacing dF. The particle dynamics is a Bloch os- 
cillation (BO) with frequency ujb — dF/fr. It should be 
mentioned that the character of these oscillations cru- 
cially depends on the type of initial conditions. We shall 
restrict ourselves by considering two limiting cases, the 
case of a completely coherent Gaussian wave-packet of 
width (To 3> 1, 

and the case of a completely incoherent wave-packet 

p(t = 0) = ^2 Pl \l){l\. (5) 

i 

In the case of coherent initial conditions ((4]) the system 
dynamics is the normal BO, where the packet center of 
gravity performs a periodic oscillation with an amplitude 
given by the Stark localization length C = J/dF, 

x(t) = C[l - cos(uj B t)} , uj B = dF/h , C = J/dF . 

(6) 

During the normal BO the wave-packet width <r(t) 
slightly oscillates, an effect which can be noticed only 
for a very weak force. This observation helps to under- 
stand the limit F — > 0, where the wave-packet spreads 
ballistically with a width a(t) obeying the equation 

cr(t) = ^Jal + (Jt/2ha ) 2 ► Jt/2ha Q , a > 1 . 

t— too 

(7) 

Note that the rate of ballistic spreading is inversely pro- 
portional to the initial width do- 
In the case of incoherent initial conditions ([5]) the Bloch 
dynamics corresponds to a so-called breathing mode, 
where x(t) = and the wave-packet width oscillates as 

o-(t) = sJal + 2C 2 sm 2 (u] B t/2) . (8) 
The F — > limit of this equation gives 

o-(i) = \Jal + 2(Jt/2h) 2 -> Jt/V2h , (9) 

where the long time limit is independent of <t . In what 
follows we shall refer to Eq. ([7]) and Eq. © as slow and 
fast ballistic regimes, respectively. 



B. Interacting atoms 

To simulate the system dynamics for a finite atom- 
atom interactions, we solve the discrete nonlinear 
Schrodinger equation, 

ifrci = - — (cj + i +ci-x) + dFlci + g\ci\ 2 ci , (10) 

where q (t) is the complex amplitude of a mini BEC asso- 
ciated with Ith well of an optical lattice and g is the ID 
macroscopic interaction constant. Following the struc- 
ture of the previous subsection we consider both coher- 
ent and incoherent initial conditions. According to ((4]), 
coherent initial conditions correspond to 

q(t = O)=v^Feaq>(*0i), (11) 

with all 9 1 — 0. To simulate the dynamics for incoherent 
initial conditions ([5]) we choose the initial phases 6i at 
random and average the solution of (flO]) over different 
realizations of random phases 6i. Typically one needs 
10 realizations to get convergence for an integrated char- 
acteristic like the wave-packet second momentum a 2 (t) 
and 100 realization to get convergence for the distribu- 
tion function Pi(t) = |cf (t)|. 

A typical weak-field evolution of the atomic density 
Pi(t) in the case of coherent initial conditions is shown 
in the upper panel in Fig. [U where time is measured in 
units of the tunneling period Tj = 2n/J. One can dis- 
tinguish several stages: The initial short-time dynamics 
corresponds to the single-particle BO, where a packet of 
interacting atoms follows the trajectory ©. This single- 
particle regime changes to the regime of dynamical in- 
stability (see, e.g., [7[) at t « Is/4, when the mean 
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FIG. 1: Subdiffusive dynamics of a BEC of interacting atoms. 
Upper panel: Evolution of the atomic density for coherent 
initial conditions, ao = 10. Lower panel: a(t) for coherent 
(solid line) and incoherent (dashed line) initial conditions. 
Parameters are J = 1, dF — 0.04, and interaction constant 
g = 10. Time is measured in units of Tj = 2n/J. 
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quasimomentum crosses the first quarter of the Brillouin 
zone. As a result the wave packet become scrambled, 
- a process which can be viewed as a formation of un- 
stable bright solitons, colliding with each other. During 
the next stage these solitons 'get thermalized' and the 
time-evolution of the site populations Pi(t) becomes a 
random process. These chaotic oscillations of the num- 
ber of atoms in any given well is a precondition for the 
subdiffusive spreading of the atomic cloud predicted in 
Ref. [13 ■ 

Before proceeding further we would like to comment 
on the relation between the mean-field treatment of the 
system, used throughout the paper, and the microscopic 
description based on the many-body Hamiltonian. With 
respect to BO this problem was addressed in the re- 
cent paper Q ■ An important conclusion one draws from 
these studies is that the discussed subdiffusive dynam- 
ics corresponds to an incoherent evolution of the single- 
particle density matrix. (Typically coherence of an initial 
BEC state is completely lost after the first 1-3 Bloch cy- 
cles.) Thus, when addressing the problem of subdiffusive 
spreading, one can use an incoherent initial state from 
the very beginning. In other words, the type of initial 
conditions affects only the transient short-time dynam- 
ics, while the long-time asymptotic dynamics is universal 
(see lower panel in Fig. [lj . Because the case of a com- 
pletely incoherent initial packet ([5]) has certain advan- 
tages from the theoretical and numerical points of view, 
in what follows we shall mainly use these initial condi- 
tions. 

Let us discuss the characteristic features of the subd- 
iffusive dynamics of interacting atoms in dependence of 
the system parameters. The dashed and solid lines in the 
lower panel of Fig. [5] show the behavior of the quantity 
@ for ,g = and g = 10, respectively. It is seen that the 
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FIG. 3: The second momentum a 2 (t) for the fixed force 
magnitude dF = 0.05 and different interaction constant 
g = 0,10,20,30,40 (from bottom to top), J = 1. The in- 
set shows o~ 2 (t) in double logarithmic scale for g = 10 and 
9 = 40. 



initial jump in the wave-packet width is due to the single- 
particle dynamics, where the maximal packet spreading 
scales as 1/F. After this jump the system enters the 
asymptotic regime, where a(t) ~ t v l 2 with v < 1. The 
upper panels in Fig. [2] depict the initial and final distri- 
bution of the site populations Pi(t). Note that this dis- 
tribution has a well-defined width and, hence, the square 
root of the second momentum is a good quantity to char- 
acterize the diffusion process [24j . Finally, Fig. [3] shows 
the behavior of the second momentum a 2 (t) for fixed F 
and different values of the interaction constant g. It is 
seen that it grows asymptotically as 
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FIG. 2: Lower panel: a(t) for g — 10 (solid line) and g = 
(dashed line). Upper panels: Initial (dashed line) and final 
(solid line) occupation probabilities in linear and logarithmic 
scale for g — 10. Incoherent initial state, J = 1 and dF = 
0.04. 



a 2 (t) 



(12) 



with an exponent v depending on g. For a large g, where 
the dynamics of the site populations is fully chaotic, the 
increment v approaches the predicted value v = 1/2 [r"^ |. 



C. Unbiased lattices 

For the sake of completeness, this subsection discusses 
the case F = 0, which requires a separate consideration. 
Indeed, for F ^ and g ^ the mean-field dynamics of 
the system is chaotic [25|]. Then, after 1-3 Bloch cycles, 
the system forgets about its initial state and enters the 
universal asymptotic regime of subdiffusive wave-packet 
spreading. However, for a vanishing static field the mean- 
field dynamics is regular and, hence, sensitive to the ini- 
tial conditions. 

It is found that in the case of a wide coherent wave- 
packet the repulsive interaction enhances the interaction- 
free slow expansion ([7]). We also mention that the 
time-evolution of the distribution function appears to be 
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FIG. 4: Ballistic spreading of interacting atoms. Lower panel: 
a(t) for F = 0, incoherent initial conditions. Values of the in- 
teraction constant g = 0, 10, 20, 30, 40 from top to bottom. 
Upper panels: Occupation probabilities at t = 207r for g = 
(left) and g — 40 (right), as compared to the initial distribu- 
tion. 



rather sensitive to the particular shape of the initial wave- 
packet. For example, for a Thomas-Fermi initial profile, 
the evolution of Pi(t) differs essentially from that for a 
Gaussian initial profile even if the wave-packet widths are 
the same. 

Unlike the case of coherent initial conditions, a repul- 
sive interaction suppresses the interaction-free fast bal- 
listic regime © [see Fig. 0] . For moderate values of the 
interaction constant, g < 10, the characteristic shape and 
time-evolution of the distribution function Pi(t) resem- 
ble those for g = 0, which is shown in the upper-left 
panel in Fig. 2J However, for stronger interactions one 
observes a qualitative deviation from the depicted shape, 
- the distribution function develops a peak at the origin 
(compare the upper-right panel) . This is a manifestation 
of the well-known phenomenon of self-trapping, where 
the system forms a soliton-like state with an energy out- 
side the Bloch band (see, for example, [26[ and references 
therein). Since the atoms belonging to the soliton state 
are permanently or temporally (soliton states with a fi- 
nite live-time) exempt from ballistic spreading, this leads 
to a decrease in a(t). In what follows, we shall take into 
account the reduced rate for ballistic spreading of in- 
teracting atoms by introducing a suppression coefficient, 
C < 1, into Eq. ©: 



a(t;g^0) = C(g/J,Jt)a(t;g = 0) 



(13) 



Althou gh t he exact analytical form of C(g/J,Jt) is un- 
known [27[, it is easy to argue that it approaches zero 
when g is increased or J is decreased. We shall discuss the 
suppression coefficient in some more details in Sec. IIII CI 
devoted to ballistic spreadings of atoms in driven lattices. 



III. DRIVEN LATTICES 
A. Single-atom dynamics 

We proceed with driven lattices. Again, our particu- 
lar interest will be the cases of completely coherent and 
completely incoherent Gaussian wave-packet. In the for- 
mer case of coherent packets a generalization of Eq. © 
for the wave-packet center of mass reads [22j 



x(t) =2\ X (t)\smct>(t) 



(14) 



where |%(t) | and <f)(t) are the absolute value and the phase 
of the following complex function 

oo 

x(t) = J £ ^(^)^expH^)sin(^), 



(15) 



In the latter case of incoherent wave-packets we have 
x(t) = and the wave-packet width oscillates as 



(16) 



A characteristic feature of the displayed equations are 
resonances at integer values of the ratio of the Bloch fre- 
quency ojb — dF/h to the driving frequency w [Hj]. To 
illustrate such a resonant dynamics of the system, we de- 
pict in Fig. [5] the width a of the incoherent packet at 
a finite time t — 2007T as a function of the driving fre- 
quency. In course of time the hight of each peak grows 
linearly with t, 



V2h ' 



(17) 
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FIG. 5: Wave-packet width a(t) at t — 200n as a function 
of driving frequency according to Eq. (|16p . Parameters are 
J = 2, dF = 0.5, Fu, = 1.21F, and a = 10. 
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while the peak tails show faster and faster oscillations 
with respect to w, with the envelope function approach- 
ing 



JJnjdFu/hu 
ft I Aw J 



(18) 



Although Eqs. (|14M16[) were obtained in a rather formal 
way, the physics behind these equations is quite simple. 
To gain a better insight into the near-resonant dynamics 
we shall use an approach involving the rotating-wave ap- 
proximation. For simplicity we shall restrict ourselves to 
the case from now on. 

It is convenient to present the Hamiltonian ([T]) in the 
basis of Wannier-Stark states ([3]) . We have 

H = dF'y^^ \m)m(m\ + dF^ cos(w£) \m)V mtm > (m'\ , 

m m,m' 

where V m ,m' — ( m \ (Sz 10 ' CI) l m ') are transition matrix 
elements between different levels of the Wannier-Stark 
ladder. Due to properties of the Bessel function these 
matrix elements differ from zero only if m! — m or ml — 
mil, 



Vrr, 



mS„ 



(z/2)8.. 



m',m±l 



z = J/dF . (20) 



Then, assuming F u <C F and using the rotating-wave ap- 
proximation, the quasienergy spectrum of the system is 
given by the Hamiltonian H — (dF — tku) \m)m(m\ + 
^dF u z^2 m (\m + l)(m| + h.c). Finally, introducing an 

effective static field dF — Ti(ujb — w) and effective tun- 
neling coefficient J = J(F LJ /2F) this Hamiltonian takes 
the form of a Hamiltonian for a stationary lattice, 

H = — ^(|m + l){m\ + h.c.) + dF^2 \m)m{m\ , (21) 

m m 

J=IE^ dF = hAuj. (22) 
2 F y ' 

Thus a near resonant driving couples Wannier-Stark 
states into new 'super' Wannier-Stark states with local- 
ization length 



C^J/dF ~ l/|Aw| 



(23) 



Using this analogy we conclude that a coherent wave- 
packet in a driven lattice will perform a super BO with 
a frequency Aw and an amplitude given in Eq. (|23l) . Of 
course, one gets the same result directly from Eqs. (| 144 - 
I15|) by keeping in the sum (I15|) only a single term with 
n = 1. Moreover, a comparison with these exact ex- 
pressions indicates that the next after the rotating-wave 
approximation corresponds to redefinition of J as 



J=JJi(F LU /F) 



(24) 



Thus the amplitude of the super BO is a nonlinear func- 
tion of the driving amplitude, - a phenomenon similar to 
the phenomenon of band collapse. 



200 



-200 



200 



-200 



25 



50 



75 



100 




FIG. 6: Super Bloch oscillations. Numerical simulation of 
the system dynamics for coherent (upper panel) or incoherent 
(lower panel) initial states. Parameters are J = 2, dF — 0.5, 
F u = 1.21F (hence J = 1), and hAuj = 0.02. 



As an illustration of the above analysis Fig.[5]shows the 
wave-packet dynamics for J = 2, dF = 0.5, F^ = 1.21F, 
and off-resonant driving Aw = 0.02, which correspond 
to J = 1 and dF — 0.02. The packet is seen to os- 
cillate with the usual Bloch frequency wb = dF/h and 
simultaneously show a slow BO with frequency Aw and 
essentially larger amplitude, given in Eq. ([23]) . We note 
in passing that such a super BO, depicted in the upper 
panel, has recently been observed for a BEC of cesium 
atoms in Ref. [2fj|. The breathing mode of super BO, 
depicted in the lower panel, was observed with not con- 
densed strontium atoms in the experiment [30j ] . 



B. Interacting atoms 

For off-resonant driving, the effect of atom-atom in- 
teractions on atomic dynamics was found to be equiva- 
lent to that in a stationary lattice, providing the former 
lattice is discussed in terms of its effective Hamiltonian 
(f2"Tj) . Thus, similar to the case of stationary lattices, in- 
teractions destroy the super BO after a few super periods 
T = ft/ 1 Aw | . The foremost consequence of the resulting 
incoherent dynamics is the formation of a smooth reso- 
nance peak (see Fig. [7]). Moreover, in course of time the 
peak shape starts to deviate from (fT5|) due to a slow in- 
crease in the wave-packet width (see Fig. [8]) . This change 
in the shape of the resonance peak may serve as an indi- 
cator of subdiffusive dynamics. 

It is interesting to study the subdiffusion with re- 
spect to the phenomenon of band collapse, described by 
Eq. (|24p. For this purpose we simulate the system dy- 
namics for finite detuning Aw, finite interaction constant 
<?, and different driving amplitudes F u . It is seen in Fig. [3] 
that at zeros of the Bessel function, where J = 0. the 
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FIG. 7: Formation of the resonance peak at ui — uib in the 
case of interacting atoms, g = 40. The other parameters are 
as in Fig. [6] (incoherent initial conditions). 
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FIG. 8: The shape of the resonance peak at t — 2007T. The 
smooth red line corresponds to Eq. (|16[) and the two broken 
lines show results of direct numerical simulations for g — 10 
(blue) and g = 40 (magenta). 

subdiffusive spreading is suppressed almost completely. 



C. Ballistic regime 

Finally we discuss the case of an exact resonance. For 
u = lob the AC field couples localized Wannier-Stark 
states into extended states and the wave-packet spread- 
ing is ballistic. However, according to results of Sec. Ill C| 
this ballistic regime may be different for different types 
of initial conditions. We have found that the case of inco- 
herent initial conditions (random phases of the complex 
amplitudes) is well captured by the effective model (p?Tj) 
but the effective model fails to describe the dynamics of 



200n 




FIG. 9: Time evolution of the width a{t) for g — 40 and 
different F w . The other parameters are J = 2, dF — 0.5, and 
HAw = 0.02. 

the original system for a coherent initial state, which is 
actually realized in a laboratory experiment. The reason 
for this is that the effective model focuses on super BO 
and overlooks ordinary BO, which appears to be impor- 
tant specifically at exact resonance. In what follows we 
analyze this situation in some more detail. 

Let us assume for the moment a non-interacting case. 
For g — and uj — lub the coherent wave packet performs 
a normal BO with slowly increasing packet width, which 
is well approximated by Eq. ([7]) for slow ballistic spread- 
ing (solid red line in the upper panel in Fig. I10[) . How- 
ever, the ordinary BO may be also dynamically unstable 
if g 7^ 0. This dynamical instability leads to an expo- 
nentially fast randomization of the relative phases of the 
complex amplitudes ci(t). As soon as the phases become 
randomized, the slow ballistic regime ([7]) changes to the 
fast one (see the upper panel in Fig. [10]) - Thus the fast 
ballistic spreading seems to be a generic case in driven 
lattices, independent of the type of initial conditions [3lj |. 

Using results of Sec. Ill C| the fast ballistic spreading is 
given by Eq. (fT3|) . where one should substitute the hop- 
ping matrix element J by the effective hopping matrix 
element (|24)l . Note that for a small J the suppression 
coefficient C ~ J and, hence, close to the zeros of the 
Bessel function the spreading rate scales as 

which should be opposed to the scaling da/dt ~ 
\Ji(Fu/F)\ for vanishing interactions. This effect is il- 
lustrated in the lower panel in Fig. [TQl showing the width 
a(t) at finite time t = 2007T in dependence on the mag- 
nitude of the driving force. The two broken curves in 
the figure are results of numerical simulations of the sys- 
tem dynamics for g = 10 and = 40, and the dashed red 
curve reproduces the dependence ^ with J substituted 
by J. A qualitative agreement with experimental results 
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as a function of time for J = 1, incoherent initial condi- 
tions, and different values of the interaction constant g. 
It is seen that for short times the spreading of interact- 
ing atoms is essentially suppressed as compared to the 
interaction-free case. Thus one may speak about tempo- 
ral self-trapping. For long times this effect of interactions 
vanishes and the expansion regime becomes ballistic, i.e., 
a(t) ~ t. These conclusions are consistent with the result 
depicted in the right- upper panel in Fig. [4j showing char- 
acteristic density profile for g ^ 0. Loosely speaking this 
density profile resembles a melting piece of ice, where the 
spreading becomes ballistic when ice melts completely. 



IV. CONCLUSION 



FIG. 10: Ballistic spreading for Auu — 0. Upper panel: The 
width o(t) for g = (red line), g = 10 (blue), and g = 
40 (magenta) and coherent initial conditions. The dashed 
red line indicates the interaction-free fast ballistic regime for 
incoherent initial conditions. Lower panel: The width a it) at 
t = 2007T. The broken lines show results of direct numerical 
simulations for g = 10 (blue), and g — 40 (magenta) for 
incoherent initial conditions. The dashed red line corresponds 
to Eq. with J substituted by J. 
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FIG. 11: Suppression coefficient C(t) = cr(t;g ^ 0)/a(t;g = 
0) as functions of time for J = 2, dF — 0.5, F u = 1.21F, and 
interaction constant g = 5, 10, 20, 30, 40 (from top to bottom). 
Incoherent initial conditions. 



In conclusion, we have studied the dynamics of inter- 
acting cold atoms in a driven tilted optical lattice (JTJ. 
Using the mean-field approach of the discrete nonlinear 
Schrodinger equation we were able to reproduce results of 
laboratory experiments [l3l |20| , which report a resonant 
response of the system at driving frequencies w ss Ws/n, 
where ljb — h/dF is the Bloch frequency. 

Our contribution in understanding the experimental 
findings is as follows. It is shown that the atomic dy- 
namics in the vicinity of any of the resonances (in this 
paper we focused on the main resonance lj w lob) can 
be described in terms of the effective Hamiltonian (|21l) . 
which formally corresponds to a stationary lattice with 
an effective static force dF = Tl(uj — ujb) and renormalized 
hopping matrix element J = JJiiF^/F). Additionally, 
one has to substitute the initial coherent state, which 
corresponds to a BEC of atoms in a laboratory experi- 
ment, by an incoherent state. The physics behind this 
seemingly artificial change of initial conditions is that, 
in a driven lattice, a BEC state rapidly decoheres in the 
presence of atom-atom interactions. 

After reformulation of the problem in terms of an ef- 
fective Hamiltonian, the atomic dynamics in a driven 
lattice can be mapped to that in a stationary lattice. 
Namely, for off-resonant driving one observes a subdiffu- 
sive spreading of the atomic cloud, where the cloud width 
grows ~ i 1 / 4 . For resonant driving one meets the sup- 
pressed ballistic spreading (|25|) , where the cloud width 
grows ~ t. 



[13|, l21| is noticed. 

For a quantitative comparison with experimental data 
a more detailed analysis of the interaction-induced sup- 
pression of ballistic spreading is needed. Figure [TT] shows 
the suppression coefficient C(t) = a(t; g ^ 0)/a(t; g = 0) 
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